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ABSTRACT 

The major theoretical hmitation for extracting cosmological parameters from the CMB 
sky lies in the precision with which we can calculate the cosmological recombination 
process. Uncertainty in the details of hydrogen and helium recombination could ef- 
fectively increase the errors or bias the values of the cosmological parameters derived 
from the Planck satellite, for example. Here we modify the cosmological recombination 
code Recfast by introducing one more parameter to reproduce the recent numeri- 
cal results for the speed-up of the helium recombination. Together with the existing 
hydrogen fudge factor, we vary these two parameters to account for the remaining 
dominant uncertainties in cosmological recombination. By using the COSMOMC code 
with Planck forecast data, we find that we need to determine the parameters to better 
than ten per cent for Hei and one per cent for H, in order to obtain negligible effects 
on the cosmological parameters. For helium recombination, if the existing studies have 
calculated the ionization fraction to the 0.1 per cent level by properly including the 
relevant physical processes, then we already have numerical calculations which are 
accurate enough for Planck. For hydrogen, setting the fudge factor to speed up low 
redshift recombination by 14 per cent appears to be sufhcient for Planck. However, 
more work still needs to be done to carry out comprehensive numerical calculations 
of all the relevant effects for hydrogen, as well as to check for effects which couple 
hydrogen and helium recombinaton through the radiation field. 

Key words: atomic processes - cosmology: cosmic microwave background ~ cosmol- 
ogy: cosmological parameters - cosmology: early universe - cosmology: observations 



1 INTRODUCTION 



Planck l|The Planck Collaborationll2006f ). the third genera- 
tion Cosmic Microwave Background (CMB) satellite will be 
launched in 2008; it will measure the CMB temperature and 
polarization anisotropies Ce at multipoles £ = 1 to ~ 2500 
at much higher precision than has been possible before. In 
order to interpret these high fidelity experimental data, we 
need to have a correspondingly high precision theory. Under- 
standing precise details of the recombination history is the 
major limiting factor in calculating the Ce to better than 
1 per cent accuracy. An assessment of the level of this un- 
certainty, in the context of the expected Planck capabilities, 
will be the subject of this paper. 

The general physic al picture of co smolo gical recombina- 
tion w as first given by iPeebled l|l968l ) and IZel'dovich et al.l 
l| 19681 ). They adopted a simple three- level atom model for 
hydrogen (H), with a consideration of the Ly a and lowest or- 
der 2s-ls two-photon rates. Thirty years later. ISeager et al.l 
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([20001) performed a detailed calculation by following all the 
resonant transitions and the lowest two-photon transition 
in multi-level atoms for both hy drogen and h elium in a 
blackbody radiation background. iLewis et al.l (|2006l ) first 
discussed how the uncertainties in recombination might bias 
the constraints on cosmological parameters coming from 
Planck; this study was mainly motivated by the efi'ect of in- 
cluding the semi-forbidden and high- order two-photon tran- 
sitions (|Dubrovich fc Grachev! |2005| ). which had been ig- 
nored in earlier calculations. 

There have been many updates and improvements in 
the mo delling of recombination since then. lSwitzer fc Hiratal 
(|2007al ) presented a multi-level calculation for neutral he- 
lium (He l) recombination including evolution of the radia- 
tion field, which had been assumed to be a perfect black- 
body in previous studies. Other issues discussed recently 
include the c ontinuum opacity due to neutral hydrogen 
(H) (see also iKholupenko et al.ll2007^ . the semi-forbidden 
transition 2'^p-l^ s (the possible importance of which was 
first proposed bv lDubrovich fc Grachevll2005l ). the feedback 
from spectral distortions between 2^p-l^s and 2^p-l^s lines, 
and the radiative line transfer. In particular, continuum 
absorption of the 2^p-l^s line photons by neutral hydro- 
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gen (H l) causes helium recombin ation to end earlier than 
previously estimated (see Fig. [1]). iHirata fc Switzeij (|2007|) 
also found that the high order two-photon rates have a neg- 
ligible effect on He I, and the same c onclusion was made 
by other groups for hydro gen as well (|Wong fc Scotti 120071 : 
IChluba fc Sunvaey l2007br). largely becaus e the a pproximate 
rates adopted by iDubrovich fc GrachevI (|2005|) had been 
overestimated. The biggest remaining uncertainty in He I 
recombination is the rate of the 2"^p-l^s transition, which 
causes a variatio n equal to about 0.1 per cent in the ioniza- 
tion fraction Xc fSwitz er fc Hiratal 2007bl ). 

For hydrogen, Chl uba et al.l ( 20071 ) improved the multi- 
level calculation by considering seperate angular momentum 
i states. This brings about a 0.6 per cent change in x^ at the 
peak of the visibility function, and about 1 per cent at red- 
shifts z < 900. The effect of the induced 2s-ls two-photon 
rate due to the radiation background (iChluba fc Su nvaev 
I2OO61 ) is p artially compensated by the f eedback of the Lya 
photons iJKholupenko fc IvanchiklbOOd ). and the net max- 
imum effect on x,, is only 0.55 per cent sd z ~ 900. 
The high-order two photon transitions^ bring about a 0.4 
per cent change in Xe at z ~ 1160 (_Wpng fc Sc otti 120071 : 
IChluba fc SunvaevI l2007bl ). There are also 0.22 per cent 
changes in ic at z ~ 1050 when one considers the Lyman 
series feedback up to n = 30, and there is additionally pos- 
sibility of direct recombinat ion, although this has on ly a 
roughly 10"* per cent effect JChluba fc SunvaevllJOOTal l. 

The list of suggested updates on Xe is certainly not com- 
plete yet, since some additional effects, such as the conver- 
gence of including higher excited states and the feedback- 
induced corrections due to the He I spectral distortions, may 
enhance or cancel other effects. In general we still need to 
develop a complete multi-level code for hydrogen, with de- 
tailed interactions between the atoms and the radiation field. 
However, what is really important here is establishing how 
these effects propagate into possible systematic uncertain- 
ties in the estimation of cosmological parameters. 

Since the un certainties in cosni ological recombination 
discussed in the iLewis et al.l (|200y) paper have been re- 
duced or updated, it is time to revisit the topic on how 
the new effects or remaining uncertainties might affect the 
constraints on cosmological parameters in future experi- 
ments. T he recent full version of t he He I recombination cal - 
culation (jSwitzer fc HiratallJOOTal lbl: iHirata fc Switzeijl2007l l 
takes too long to run to be included within the current 
Boltzmann codes for Ci. So instead, in this paper, we try 
to reprod uce the updated i onization history by modifying 
Recfast jSeager et al.l 19991 ) using a sim ple parametrization 
based on the fitting formulae provided by_'KholuDcnko ct al. 
|2003). We then use the COSMOMC f|Lewis fc Bridle 2002 ) 
code to investigate how much this impacts the constraints 
on cosmological parameters for an experiment like Planck. 
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Figure 1. Top panel: Ionization fraction Xc as a function of 
redshift z. The dotted (red) line is calculated using the origi- 
nal REC FAST code. The solid (bl ack) line is the numerical re- 
sult from lSwitzer fc Hiratal pOOTo ). while the dashed (blue) and 
long-dashed (gree n) lines are both e valuated based on the modi- 
fication given by iKholupenko et al.l 112007;) - the dashed one has 
''He = 0.97 (the value used in the original paper) and the long- 
dashed one has ftne = 0.86. Bottom panel: The visibility func- 
tion g{z) versus redshift z. The two curves calculated (dotted 
and long-dashed) correspond to the same recombination models 
in the upper panel. The cosmological parameters used for these 
two graphs are, Ob = 0.04592, Hm = 0.27, OcDM = 0.22408, 
TcMB = 2.728 K, Ho = 71 kms"! Mpc"! and /hc = 0.079. 
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Figure 2. Ionization fraction Xc as a function of redshift z cal- 
culated based on the modified He I recombination discussed here 
with different values of the helium fitting parameter feno- The 
curve with 6hc = (dot-dashed, blue) overlaps the line using 
Saha equilibrium recombination (solid, black). The cosmological 
parameters used in this graph are the same as for Fig. [T] 



2 RECOMBINATION MODEL 

In this paper, we modify Recfast based on the fitting for- 
mulae given by iKholupenko et al.l (|2007l ) for including the 
effect of the continuum opacity of neutral hydrogen for He I 
recombination. The basis set of rate equations of the ioniza- 



tion fraction of H and He I used in Recfast are: 

i/(.)(l + .)5i = 
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the He I atom with electron in the 2 s state, and /ii'23a,c is 
the ionization energy of the 2^s state. 

The correction factor ChgI accounts for the slow recom- 
bination due to the bottleneck of the He I 2^p-l^s transi- 
tion among singlets. We can also derive the corresponding 
correction factor Chci for the triplets. The Ku, Kuei and 
-R'hoi quantities are the cosmological redshifting of the H 
Lya, He I 2^p-l^s and He I 2''p-l^s transition line photons, 
respectively. The factor K used in Recfast is a good ap- 
proximation when the line is optically thick (r 3> 1) and 
the Sobolev escape probability ps is roughly equal to 1/r. 
In general, we can relate K and ps through the following 
equations (taking He I as an example): 



Figure 3. The ionization fraction Xc as a function of redshift z 
calculated with different values of the hydrogen fudge factor Fjj . 
The cosmological parameters used in this graph are the same as 
in Fig. [1] 
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Note that x^ is defined as the ratio of free electons per H 
atom and so a;c > 1 during He recom bination. We follow the 
exact notation used in iSeager et al.l (|.1999) and we do not 
repeat the definitions of all symbols, except those that did 
not appear in that paper. The last term in equation ([2]) is 
added to the original dxHeii/dz rate for the recombination 
of He I through the triplets by including the semi-forbidden 
transition from the 2^p state to the f^s ground state. This 
additional term can be easily derived by considering an extra 
path for electrons to cascade down in He I by going from 
the continuum through 2^p to ground state, and assuming 
that the rate of change of the population of the 2^p state 
is negligibly small. The superscript 't' stands for triplets, 
so that, for example, Qnei is the Case B Hei recombination 
for tri plets. Based on the data given bv lHummer fc StorevI 
(|l99g ). ajiei is fitted with the same fun ctional form used fo r 
the and singlets (see equation (4), in ISeager et al.lll999r ). 
with different values for the parameters: p = 0.761; q = 
^Q-16.306. y^ ^ ;^q5.ii4 j^. ^^^ Ta = 3K. This fit is accurate 

to better than 1 per cent for temperatures between 10^'^ 
and 10* K. Here Pli^i is the photoionization coefficient for 
the triplets and is calculated from QhcI by 
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where Phc+ ^^'^ 3hcI,23s sxe the degeneracies of He"*" and of 
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where ^Hci,2ip-iis ^^'^ ^Hci,23p-iis s-re the Einstein A co- 
efficients of the He I 2'^p-l^s and He I 2^p-l'^s transitions, 
respectively. Note that Afj^j y3p_Ti „ = gHpr . c^^p^/flHr r 



A 



HcI,23Pi-lls 



= 1/3 X 



177.58 8^1 l|Lach fc PachuckilboOll ) 



For He I 2^p-l^s, we replace ps by the new escape prob- 
(2) ability Posc, to include the effect of the continuum opac- 
ity diw_to_Jfj_based_on_tlie approximate formula suggested 
bv lKholupenko et aj] (|2007l '). Explicitly this is 



Pesc = PS + Pcon.H , 



where 



PS = 



Pcon.H ~ 



and 



1 



with 



1 + auoT^Ho 



ShoI.iIs ^HcI 
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Ho — a;Hcii)c 



87r3/2crH,ls(j^HcI,2lp)^'HcI,2lp^^D.2lp(l " ^p) ' 

where crH,is(j^Hci,2ip) is the H ionization cross-section at fre- 
quency i-'HoI,2lp and Al/D_2lp = J^Hci,2ip\/2fcBrM/mHcc2 is 
the thermal width of the He I 2'^p-l^s line. The 7 factor in 
Pcon,H is approximately the ratio of the He I 2^p-l^s transi- 
tion rate to the H photoionization rate. When 7^1, the 
effect of the continuum opacity due to neutral hydrogen on 
the He I recombination is negligible. Here ahc and 6hc are 
fitting parameters, whi ch are equal to 0.36 a nd 0.97, based 
on the results from Kh olupenko et al.l (|2007l ). 

We now try to reproduce these results with our modi- 
fied Recfast. Fig. [l](upper panel) shows the numerical re- 
sult of the ionization fraction Xc fronr di fferent He I recom- 
bination calculatio ns. The re s ults fr om iKholupenko et al.l 
(2007) and Switze r~fc"Hiratal (|2007bh both demonstrate a 
significant speed up of He I recombination compared with 
the original Recfast. We do not expect th e se tw o curves 
to match each other, since Kholupenko et al.l (|2007l ) just in- 
cluded the effect of the continuum opacity due to hydro- 
gen, which is only one of the main improvements stated in 
Switzer fc Hirata ll2007bl ). Nevertheless, we can regard the 
Kholupenko et al.l (|2007l ) study as giving a simple fitting 



model in a three-level atom to account for the speed-up of 
the He I recombination. Fig.[2]shows how the ionization his- 



4 W. Y. Wong, A. Moss and D. Scott 




0.022 0.0225 0.023 0.1 0.105 



1.036 1.039 1.04 







0.06 0.06 0.1 0.12 0.93 0.94 0.95 0.96 0.97 



logllO'" A.] 







Figure 4. Marginalized posterior distributions for forecast 
Planck data varying the hydrogen recombination only. All the 
curves are generated using the original Recfast code. The 
solid (black) curve uses fixed Fjj , while the dotted (red) and 
dashed (green) allow for varying F^ with Gaussian distributions 
centred at 1.14, with a = 0.1 and 0.01, respectively. Note that us- 
ing a flat prior (between and 1.5) for Fjj gives the same spectra 
as the case with cr = 0.1 (the red dotted line). 



tory changes with different values of the fitting parameter 
bno (with auc fixed to be 0.36). When bua is larger than 
1.2, the effect of the neutral H is tiny and the fit returns to 
the situation with no continuum opacity. However, if 6hc is 
smaller than 1, the effect of the continuum opacity becomes 
more significant. Note that when 6hc = 0, both the escape 
probability Pcsc and the correction factor ChoI are close to 
unity. This means that almost all the emitted photons can 
escape to infinity and so the ionization history returns to 
Saha equilibrium for He l recombination. 

This simple fitting formula can reproduce quite well the 
detailed numerical result for the ionization history at the 
later stages of Hel recombination. From Fig. [l] we can see 
that our model wit h fene = 0-86 matches wit h the numerical 
result at z < 2000 (|Switzer fc HiratallJOOTbl ). Although our 
fitting model does not agree so well with the numerical re- 
sults for the earlier stages of He I recombination, the effect 
on the Ce is neligible. This is because the visibility function 
g{z) = e'^dr/dz, is very low at z > 2000 (at least 16 orders 
of magnitude smaller than the maximum value oi g{z)), as 
shown in the lower panel of Fig. [T] Our fitting approach also 
appears to work well for other cosmological models (Switzer 
& Hirata, private communication). 

In this paper, we employ the fudge factor Fh for 
H (which is the extra factor multiplying qh) and the Hel 
parameter 6hc in our model to represent the remaining un- 
certainties in recombination. For Hel, the factors auc and 
6He in equation (|ll|l are quite degenerate. We choose to fix 
aHe and use fenc as the free parameter in this paper; this 
is because it measures the power dependence of the ratio of 
the relevant rates 7 in the escape probability due to the con- 
tinuum opacity Pcon,H- For hydrogen recombination, all the 
individual updates suggested recently give an overall change 
less than 0.5 per cent in x^ around the peak of the visibility 
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Figure 5. Marginalized posterior distributions for forecast 
Planck data with hydrogen and helium phonomenological param- 
eters both allowed to vary. The solid (black) curve shows the con- 
straints using the original Recfast code and allowing Fn to be 
a free parameter. The other curves also allow for the variation of 
i^H and use the fitting function for Hel recombination described 
in Section 2: the dotted (green) line sets 6hc equal to 0.86; the 
dashed (red) one is with a fiat prior for bjje from to 1.5; and the 
long-dashed (blue) one is with a narrow prior for 6hci consisting 
of a Gaussian centred at 0.86 and with a = 0.1. 



function. Only the effect of considering the separate ^-states 
causes more than a 1 per cent change, and only for the final 
stages of hydrogen recombination [z < 900). Therefore, we 
think it is sufficient to represent this uncertainty with the 
usual fudge factor Fh, which basically controls the speed of 
the end of hydrogen recombination (see Fig. [3]). The best-fit 
to the current recombination calculation has i*' ~ 1.14. 



3 FORECAST DATA 



We use the COSMOMC (iLewis fc Bridle! |2002| ) code to per- 
form a Markov Chain Monte Carlo (MCMC) calculation for 
sampling the posterior distribution with given forecast data. 
The simulated Planck data and likeliho od function are gen - 
erated based on the settings suggested in lLewis et al.l l|2006l ). 
We use full polarization information for Planck by consider- 
ing the temperature T and i?-type polarization anisotropies 
for £ ^ 2400, and assume that they are statistically isotropic 
and Gaussian. The noise is also isotropic and is based on 
a simplified model with iVf^ = Nf^/i = 2 x 10-''^K^ 
having a Gaussian beam of 7 arcminutes (Full Width Half 
Maximum. lThe Planck Collaborationll2006l ) . For our fiducial 
model, we adopt the best values of the six cosmological pa- 
ram eters in a ACDM mo del from the WMAP three-year re- 
sult (|Spergel et al.ll2007l ) . The six parameters are the baryon 
density Qhh^ ~ 0.0223, the cold dark matter density Q,ch^ = 
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0.104, the present Hubble parameter Ho = 73 kms ^Mpc ^ , 
the constant scalar adiabatic spectral index Us — 0.951, the 
scalar amplitude (at k = 0.05 Mpc"^) 10^° A^ = 3.02 and the 
optical depth due to reionization (based on a sharp transi- 
tion) T = 0.09. For recombination, we calculate the ioniza- 
tion history using the original Recfast with the fiidge factor 
for hydrogen recombination Fu set to 1.14 and the helium 
abundance equal to 0.24. 

In this study, we only vary the basic six standard cosmo- 
logical parameters stated above, together with the hydrogen 
fudge factor Fu and also the helium 6ho factor for the re- 
combination process. Of course degeneracies will in general 
be worse if one allows for a wider set of parameters. 



4 RESULTS 

Fig. U] shows the parameter constraints from our forecast 
Planck likelihood function using the original Recfast code 
with varying Fh and adopting different priors. For the 
Planck forecast data, Fu can be well c onstrained away from 
zero (the same result as in lLewis et al.ii 2006 ') and is bounded 
by a nearly Gaussian distribution with a approximately 
equal to 0.1. When we only vary Fh with different pri- 
ors (compared with fixing it to 1.14), it basically does not 
change the size of the error bars on the cosmological parame- 
ters, except for the scalar adiabatic amplitude lO^^As. From 
Fig. [3] we can see that the factor Fh controls the speed of 
the final stages of H recombination, when most of the atoms 
and electrons have already recombined. Changing Fh affects 
the optical depth r from Thomson scattering, which deter- 
mines the overall normalization amplitude of Ci (oc e"^'^) 
at angular scales below that subtended by the size of the 
horizon at last scattering (£ > 100). This is the reason why 
varying Fh affects the uncertainty in As, since As also con- 
trols the overall amplitude of Ce (see the upper right panel 
in Fig. |6] for the marginalized distribution for Fh and As). 
The modified recombination model also changes the peak 
value (but not really the width) of the adiabatic spectral in- 
dex ris distribution, as one can see by comparing the dotted 
and dashed curves in Fig. 3] 

Based on all the suggested effects on H recombination, 
the uncertainty in Xc is at the level of a few per cent at 
z < 900, which corresponds to roughly a 1 per cent change 
in Fh. In Fig. 2] we have also tried to take this uncertainty 
into account by considering a prior on Fh with a = 0.01 
(the long-dashed curves). We find that the result is almost 
the same as for the case using ct = 0.1 for the Fh prior. On 
the other hand, the error bar (measured using the 68 per 
cent confidence level, say) of As is increased by 40 and 16 
per cent with a — 0.1 and 0.01, respectively. 

Fig. [5] shows the comparison of the constraints in the 
original and modified versions of Recfast, with both H and 
He parameterized. By comparing the solid and dotted curves 
in Fig. [5] we can see that only the peaks of the spectra of the 
cosmological parameters are changed, but not the width of 
the distributions, when switching between the original and 
modified Recfast codes. Allowing buc to float in the mod- 
ified recombination model only leads to an increase in the 
error bar for spectral index Wg among all the parameters, in- 
cluding Fh. For the dashed curves, we used a very conserva- 
tive prior for 6hg, namely a flat spectrum from to 1.5 (i.e. 
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Figure 6. Projected 2D likelihood for the four parameters ris, 
As, Fh and feno- Shading corresponds to the marginalized proba- 
bilities with contours at 68 per cent and 95 per cent confidence. 



from Saha recombination to the old Recfast behaviour). 
We can see that the value of 6hc is poorly constrained, be- 
cause the CMB is only weakly sensitive to the details of He I 
recombination. Nevertheless, this variation allows for faster 
He I recombination than in the original Recfast code and 
this skews the distribution of ris towards higher values (see 
also the upper left panel in Fig. IS}. This is because a faster 
He I recombination leads to fewer free electrons before H re- 
combination and this increases the diffusion length of the 
photons and baryons. This in turn decreases the damping 
scale of the acoustic oscillations at high £, which therefore 
gives a higher value of ris. In addition, this variation in 6hc 
increases the uncertainty (at the 68 per cent confldence level) 
of ris by 11 per cent. 

Ba sed on the comprehensive study of ISwitzer fc Hiratal 
(J2007b|) , the dominant remaining uncertainty in He I recom- 
bination is the 2^p-l^s transition rate, which causes about 
a 0.1 per cent variation in Xc ai z ~ 1900. For our fltting 
procedure this corresponds to about a 1 per cent change in 
bua- We try to take this uncertainty into account in our cal- 
culation by adopting a prior on &hc which is peaked at 0.86 
with width (sigma) liberally set to 0.1. From Fig.[S] one can 
see that the error bar on ris is then reduced to almost the 
same size as found when fixing 6hc equal to 0.86 (the dotted 
and long-dashed curves). This means that, for the sensitiv- 
ity expected from Planck, it is sufficient if we can determine 
&Hc to better than 10 per cent accuracy. 

As well as the individual marginalized uncertainties, we 
can also look at whether there are degeneracies among the 
parameters. From Fig. |S] we see that Fh and 6hc are quite 
independent. This is because the two parameters govern re- 
combination at very different times. As discussed before, 6hc 
controls the speed of He I recombination, which affects the 
high-z tail of the visiblity function, while Fh controls the 
low- 2: part. 
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5 DISCUSSION AND CONCLUSIONS 

In this paper, we modify Recfast by introducing one more 
parameter &hc (besides the hydrogen fudge factor _Fh) to 
mimic the recent numerical results for the speed-up of He I 
recombination. By using the COSMOMC code with forecast 
Planck data, we examine the variation of these two factors 
to account for the remaining dominant uncertainties in the 
cosmological recombination calculation. For He I, the main 
uncert ainty comes from the 2^p-l^s rate IjSwitzer fc Hiratal 
[2007bj), which corresponds to about a 1 per cent change 
in bue- We find that this level of variation has a negligible 
effect on the determination of the cosmological parameters. 
Therefore, based on this simple model, if the existing studies 
have properly considered all the relevant physical radiative 
processes in order to provide x^ to 0.1 per cent accuracy 
during He I recombination, then we already have numerical 
calculations which are accurate enough for Planck. 

For H, since there is still no comprehensive model which 
considers all the interactions between the atomic transitions 
and the radiation background, we consider the size of the 
updates as an indication of the existing level of uncertainty. 
We represent this uncertainty by varying the fudge factor 
Fh, because the largest update on Xb occurs at z < 900, and 
comes from a consideration of t he separate angular momen- 
tum states jChluba et al.l 120071 '). We find that Fn needs to 
be determined to better than 1 per cent accuracy in order to 
have negligible effect on the determination of cosmological 
parameters with Planck. 

Hydrogen recombination is of course important for the 
formation of the CMB anisotropies Ct, since it determines 
the detailed profile around the peak of the visibility function 
g{z). A comprehensive numerical calculation of the recombi- 
nation of H (similar to He l) to include at least all the recent 
suggestions for updates on the evolution of Xe is an urgent 
task. We need to determine that the phenomenological pa- 
rameters Fu and bue are fully understood at the < 1 per 
cent level before we can be confident that the uncertainties 
in the details of recombination will have no significant ef- 
fect on the determination of cosmological parameters from 
Planck. 



Chluba J., Sunyaev R. A., 2007, ,arXiv:0705.3033l 
Dubrovich V. K., Grachev S. I. 2005, Astronomy Letters, 

31, 359, arXiv:astro-ph/0501672' 

Hirata, C. M., Switzer, E. R. 2007, arXiv:astro-ph/0702144] 
Hummer D. G., & Storey P. J. 1998, MNRAS, 297, 1073 
Kholupenko, E. E., Ivanchik, A. V. 2006, Astronomy Let- 
ters, 32, 795 
Kholupenko E. E., Ivanchik A. V., Varshalovich D. A. 2007, 

MNRAS, L42 
Lach G., Pachucki K. 2001, Phys. Rev. A, 64, 042510 
Lewis A., Bridle S. 2002, Physical Review D, 66, 103511 
Lewis A., Weller J., Battye R. 2006, MNRAS, 373, 561 
Peebles P. J. E. 1968, Astrophysical Journal, 153, 1 
The Planck CoUaboration 2006, ESA-SCI(2005)1, 



arXiv:astro-ph/0604069 



Scott, D. 1999, The Astrophysical 
2000, Astrophysical 



Seager S., Sasselov D. D. 

Journal, 523, LI 
Seager S., Sasselov D. D., Scott D. 
Journal Supplement Series, 128, 407 
Spergel D. N., et al. 2007, Astrophysical Journal Supple- 
ment, 170, 377 

Switzer E. R., Hirata, C. M. 2007, |arXiv:astro-ph/0702143| 
Hirata C. M. 2007, arXiv:astro-ph/0702145 
Ferland G. J. 1996, Astrophysical Journal 
103, 467 

Scott D. 2007, MNRAS, 375, 1441 
B., Kurt V. G., Syunyaev R. A. 1968, Zh. 
Fiz., 55, 278; English translation, 1969, Soviet 



Switzer E. R., 
Verner D. A., 

Supplement, 
Wong W. Y., 
Zel'dovich Y. 

Eksp. Teor. 

Phys. -JETP Lett., 18, 146 



This paper has been typeset from a TjrjX/ LT^jX file prepared 
by the author. 



6 ACKNOWLEDGEMENTS 

This work made use of the code CoSMOMC, written by 
Antony Lewis and Sarah Bridle. We are also grateful to 
Chris Hirata and Eric Switzer for comparing our results 
with their own calculations. This research was supported by 
the Natural Sciences and Engineering Research Council of 
Canada and the Canadian Space Agency as well as through 
a University of British Columbia Graduate Fellowship. A 
new version of the Recfast code (version 1.4) is available 
from ^ http://www.astro.ubc.ca/people/scott/recfast. html. 



REFERENCES 

Chluba J., Sunyaev R. A. 2006, Astronomy and Astro- 
physics, 446, 39 

Chluba J., Rubino-Martin J. A., Sunyaev R. A. 2007, MN- 
RAS, 374, 1310 

Chluba J., Sunyaev R. A. 2007, jarXiv:astro-ph/0702531 



